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Two devices arc proposed for measuring absorption coefficients in weakly absorbing materials. 
The first device measures cylindrical samples and the second device measures flat plate or disk 
samples. This paper reports on the derivations for the steady-state and transient solutions to the heat 
diffusion equations which describe the barothermal behavior of the two proposed devices. In addition. 
Green's function -techniques are used to describe the cyclic heating and cooling of the cylinders and 
plates. 
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1. Introduction 

Optical communications, integrated optics, and high 
power lasers are developing technologies which depend 
in part upon highly transparent solids. For example, 
optical communication systems use optical fibers with 
absorption coefficients which are less than 10~ 4 cm -1 at 
operating wavelengths [l]. 1 Continuing further ad- 
vances in these technologies and determining the ab- 
sorption mechanisms in highly transparent materials 
require improved methods for measuring very small 
absorption coefficients. 

Four general techniques for measuring absorption 
coefficients exist. These techniques involve respec- 
tively calorimetry, spectrophotometry, emissivity, and 
gas pressure measurements. 

In the calorimetric technique, one measures the 
temperature rise of the solid which occurs when a frac- 
tion of the energy in a beam of radiation is absorbed. 
When the temperature rise, reflectivity, and laser 
power are known to within an accuracy of 1 percent, 
this method is capable of giving absorption coefficients 
on the order of 10 ~ 4 cm -1 within an accuracy of a few 
percent. This method then has reasonably high sen- 
sitivity, and requires laser beam powers of the order of 
tens of mW for absorption coefficients of the order of 
10~ 3 cm -1 . It also is insensitive to those scattering 
centers which do not absorb any of the radiation. How- 
ever, the sample must be contacted with thermo- 
couples to measure its temperature and the absorption 
coefficient can be determined only at those wave- 
lengths for which sufficiently powerful lasers are 
available. 



i brackets indicate the literature references at the end of this 



The second technique involves measuring the trans- 
mission loss of a beam of radiation passing through the 
solid. Variations of this technique include either a 
single beam or double beam spectrophotometer and 
may use samples of two different lengths. Double beam 
methods using samples of different lengths give prod- 
ucts of the absorption coefficients a with the sample 
thickness d as small as ad ~ 0.002 with an accuracy of 
±0.0002. The photometric technique has the advantage 
of using light sources with continuous values of wave- 
lengths present; but it has the disadvantages that 
sample preparation is important to avoid surface 
scratches and dust and that the samples should noU 
have scattering centers such as inclusions. 

In the third technique one measures the emission of 
the sample at a given temperature and compares it 
directly to a blackbody at the same temperature. When 
the product ad is very small, the emittance E (the ratio 
of thermal radiation per unit area emitted by the 
sample to that emitted by a blackbody at the same tem- 
perature) becomes E ~ ad. For values of ad less than 
0.001, the emittance measurements can be more sensi- 
tive for determining a than the photometric measure- 
ments. In addition, by this technique, one can obtain 
values of absorption coefficients over a continuous 
range of wavelengths. But exceptional temperature 
control is essential to achieve such sensitivity. 

In the fourth technique, the energy absorbed from a 
beam of radiation passing through the sample produces 
heat. The heat then diffuses from the solid to a non- 
absorbing, confined gas which is adjacent to the sam- 
ple. The heat transfer process leads to a pressure rise 
in the gas. Because this technique tends to average 
over any absorption inhomogeneities in the sample, we 
also expect that localized absorbing or scattering 
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centers could be less of a problem for this technique 
than for the calorimetric techniques in which the 
placement of thermocouples may be important. 

To develop further this last technique, Bennett and 
Forman recently proposed alternative ways to measure 
very small absorption coefficients. They considered in 
a series of papers [2-6] both barothermal and photo- 
acoustic modes of operation for devices to characterize 
weakly absorbing materials. In these papers, they 
cited only the results for the solutions to the heat 
diffusion equations which describe the behavior of 
the proposed devices. They did not derive the solutions 
because they wanted to stress in those papers the major 
conclusions for designing absorption devices with 
optimum performance. 

Figures 1 and 2, respectively, show schematics of 
their proposed devices for measuring absorption co- 
efficients of long cylindrical samples and of thin disk or 
plate-shaped samples. Each of these devices has a 
characteristic frequency which is related to the inverse 
of the time required for transferring a heat pulse from 
the solid to the gas. The barothermal mode of opera- 
tion occurs whenever the modulation frequency of the 
beam of radiation is less than the characteristic fre- 
quency of the device; and the photoacoustic mode of 
operation occurs whenever the modulation frequency 
of the beam of radiation is greater than the character- 
istic frequency of the device [7]. 

In this paper, the author gives derivations for the 
steady-state and transient solutions to the heat diffu- 



sion equations which describe the barothermal be- 
havior of the two proposed devices. The equations 
which describe the photoacoustic behavior of the 
devices are derived in (separate papers [4, 6]). 

The steady-state and transient solutions for cylin- 
drical samples are derived here, respectively, in sec- 
tion 2 and section 3; and the steady-state and transient 
solutions for disk or plate-shaped samples are de- 
rived, respectively, in section 4 and section 5. The 
Green's functions which describe the cyclic heating 
and cooling of the cylinders and plates are constructed 
in section 6. And finally, frequently occurring integrals 
involving the eigenfunctions for the equations are given 
in appendix A for cylinders and in appendix B for disks 
or plates. 

2. Steady-State Solution — Cylinder 

From [3], the steady-state temperature v(r) satisfies 
the following heat diffusion equations for the proposed 
device in figure 1: When ^ r^ r s , 



, <Pv,(r) k, dv,(r) , <?«(r) . 

k, H 1 — = 0; 

dr- r dr p,L s 



and when r,< r«£ n, 



cP-Vfi{r) kg dvg(r) 



dr 2 



+ 



dr 



= 0. 



(1) 



(2) 




1--C/2 



i4A 



FIGURE 1. Schematic of an absorption measuring device {barothermal) for long cylindrical samples. 

Insulating annuli support coaxially the long cylindrical transparent sample of radius r s inside a larger cylinder of radius r,. A confined 
nonabsorbing gas at ambient pressure p<> fills the space between the two cylinders. The walls of the outer cylinder are at a constant tem- 
perature 7V A collimated beam of radiation, such as a laser beam, propagates coaxially through the cylindrical sample. The beam of 
radiation has a power Wi and an effective radius n. A pressure transducer, which is not shown, monitors the pressure at the heat sink — 
gas interface. 
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FIGURE 2. Schematic of an absorption measuring cell (barot hernial) for plate-shaped samples. 
\ confined, nonabsorbing gas at ambient pressure/),, tills the space inside the cefl between two windows. The walls of the cell may have 

a polygonal or circular cross section The windows which enclose the ^as at each end of the cell are made from the weakly absorbing 
material under investigation. The external faces of the two windows are at a constant temperature To. \ collimated heam of radiation, 
such as a laser heam, with power W\ and effective radius ft propagates along the Z direction through the pressure cell and its windows. 
\ pressure transducer, which is not shown, measures the pressure in the gas at \z\ =£ z 2 . 



The temperature in the solid is v s (r) and the tem- 
perature in the gas is v !f (r). The power absorbed per 
unit volume in the solid is 



Qs(r) 



Q,Un0^r 
Ofor r> n, 



r 8 



(3) 



where Qi = ct 8 Wilirrf, a s is the absorption coefficient 
of the solid, W\ is the power of the beam of radiation, 
and ri is the effective radius of the collimated beam. 
The radius of the cylindrical sample is r 8 and the radius 
of the outer cylinder (heat sink) is r>. The thermal 
diffusivity k = K/pC, where K is the thermal conduc- 
tivity, p is the density, and C is the specific heat at 
constant volume. The subscript s denotes quantities 
for the solid and the subscript g denotes quantities for 
the gas. 

The boundary conditions are that 



tfr(n) =r , 



where To is the temperature of the heat sink. Con- 
tinuity of temperature across the solid-gas interface 
at r= r s gives, 

v 8 (r 8 )=v g (r g ), (5) 

and conservation of heat flow across the solid-gas 



interface at r= r 8 gives 



K s 



dr 



dv<, 



— K<i 
>•=/• dr 



(6) 



Rearranging eq (1) and eq (2) yields, respectively, 



and 



drV dr ) K, 



i{'& 



(7) 



(8) 



and integrating these equations twice yields the steady 
state solutions; 



(4) v 8 (r)=v 8 (0)- 



Q[ [r 2 forO^r^r, 

4iC*lrf[l + 21ii (r/r,)] for n^r^r 8> 



(9) 



and 



v g (r) = Cigln r-\-C 2 <, for r s ^ r^ n , (10) 



where C\ u and C%g are constants of integration and 
where the limit that 
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t • dv * n 

Lim r -7- = 

dr 



r-^0 



(ID 



and the temperature continuity equation at r=r s . 
requires that 



is invoked to obtain eq (9) (i.e., that v s remains finite at 
r=0). 

Boundary conditions (4), (5), and (6) require that the an d 
temperature at the center of the cylinder r=0 is 



R,(r s )=R a (r t ) 



(20) 



*<0)-r.+ gUl-4S.if*y + 2]n(& 



&=(kjk.)? r (21) 



4K S 



Ka 



r> 



and that the temperature in the gas is 

Qui 



(12) Also, the constant B s = because the temperature must 
be finite for all r ^ and t. 

Hence, the conservation of heat flow across the solid- 
gas interface 



v g (r) =T 



2K g ln yn 



(13) 



3. Transient Solution — Cylinder 



The transient temperature u(r, t) for times t^ is 
given in [3] by the diffusion equations 



. d 2 , , k 8 du s (r, t) du s (r, t) 
k 8 — u 8 {r,t) H — 
dr- r 



dr 



d* 



and 



a 2 ^ , k 9 dUg(r, t) dugjr, t) 
k g r- 2 u g (ri t) + — , 

or r dr dt 



(14) 



(15) 



K s du r {r, t)ldr\r=r=Kgdu g (r, t)ldr\ r ~r 8 (22) 

and the other boundary conditions eq (19) and eq (20) 
yield three homogeneous simultaneous equations for 
the three unknowns A 8 , A g , and B g . Solutions for this 
set of three equations exist if and only if the de- 
terminant, D, of their coefficients vanishes; the vanish- 
ing of the determinant yields the eigenvalues f#=f j, 
where 7=1, 2 . . .. The 7th eigenvalue £j is that 
value of f g for which the determinant vanishes for the 
7th time, i.e., fj<£/ + i; namely, det Z)(fj)=0, where 

MUt) 

Joi^n) -Mfr.) 



D 



At time ( = 0, the sum of the transient and steady- 
state temperatures is everywhere equal to the ambient 
temperature 7\>, 



-Y»U„r s ) 

KgigYtdgT,) 



(23) 



T(r, 0) = u(r, 0)+v(r) = T ; 



(16) 



and at any time t the temperature at the heat sink 
r= n is To, i.e., 



because 



Ug(n, t) = 0, 



Lim u(r, t) = 0, 



(17) 



The eigenvalues £/ also determine the character- 
istic times t; for the solid-gas system; namely, 

tj=(£]kg)-K (24) 

Expressing^., and B u in terms oiAj — A,,, we obtain, 



/?,(r) = 2^F () fe, r u r § ){Jo(ijr)lM{jr.)h (25) 



for0 = 



by the definition of transient solutions and because 

A general separable solution to an equation having 
the form of either eq (14) or (15) is 



r ^ r, s and 

Rg(r) = ^AjFMi,rur) 



j=i 



u(r, t) =R(r)r(t) = [AM£r) + £F (£r)] 

exp {—} 2 kt). 



(18) 



for r s ^, r ^ n> where ^ l — (kjk s ) £ 2 and where the 
function F n is 

X? i tL \ 

Aylkln)] nU,r) -J^,rdYn^,r)}IY^ uri ). (27) 

Relations (25), (26), and (27) are valid only if 
kind [8], f is the constant of separation, and A and B Jo(^jr s ) # and if Yo(^jrt) ^ 0. 

are constants. The quantities A, B, £, and k have the In terms of the functions F H (£, n, r), the eigenvalue 

subscript s for the solid and the subscript gfor the gas. equation satisfied by the eigenvalue fj becomes, 
Equation (17) states that 

Krf//i(for.,.F'o(&, r u r.) -K^jJo({jr g )Fi(( , n, r.) 
« ff (r0=0, (19) =0. 



where /,, is the nth order Bessel function of the first 
kind, Y n is the nth order Bessel function of the second 
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We determine the coefficients Aj by using the t = (A4) and the fact that dBo(^r) ldr= — B i(fr) , we find 

that A(£j, £/) vanishes unless fj=£i; and when &=f< 
the factor A (fj, £/) = A ; becomes 



relations obtained from eq (16); namely 

R 8 (r) = To-v 8 (r) for O^r^r, (28) 

and 

f? // (r) = 7 1 ( ,-2;,(r) for r, ^ r ^ r*. (29) 

First, we consider the eigenfunctions for the solid- 
gas system 



£<)(£;, r« f r) 



-{ 



Fo(ejJ-i.r.){Jo(€jr)IM€jf,)} 
Fo(€j, ru r) 



for s= ; r =£ r, 

for r. v ^ r =£ r, i. 

(30) 



Next, we multiply both sides of eq (28) and eq (29) by 
K, 



and by 



Kg 



rEo(£i, ru r) whenO"^ r ^ r 8 



rEo{£u ru r) when r 8 ^ r ^ n 



and then we integrate the resulting expressions over 
the interval O^r^r,. This procedure yields the fol- 
lowing expression for the coefficients Aj 



= ^\ {r -».(0)} f "£„(£, ,r,.r)rdr 



(31) 



+^j "*.<£i, n,r)rdr 



Jri 



1 + 2 In (r/r,)}£o(fi, r,, r)r ! dr 



+ff 






if" In 



(r/r/)£ (f/, H, r)rdr 



/here 



A(&, fi 



K, 



[['EAfaru 



r)E 9 (^i, r t , r)rdr 



+ 



* f "*•(&, r«, 

g J r s 



2L 



KJi„ 
K„k s ' 



-1 \Fl(i J ,r i , r.) 



+ (»f/f*)F«(6, r f , n) + 






1)F?(&, r,, r.) 



(33) 



Applying eqs (Al) to (A3) to the right-hand side of 
eq (31), we obtain, after several steps, an expression 
for the coefficients Aj, namely, 



4 = - 



Qiri Ji(£r,)F<,(tj,ri,r,) 



(34) 



Because we sum in the computer program used in 
[3] a finite number of terms in such expressions as 
eq (25) and eq (26), we normalize them's for 1 ^y =s /V 
by using the relation (29) evaluated at r= r 8 . To accom- 
plish this normalization, we introduce the normaliza- 
tion factor a.v and replace Aj with a.\Aj{N) in the 
relation (29). 'Phis then gives us an expression for the 
factor a. v ; namely, 

aN= ~2K a ln {n )/% A lWF*{£ h r u r 9 ). (35) 



4. Steady-State Solution — Plate 

From [5], the steady-state temperature v(z) 
satisfies the following heat diffusion equations for 
the proposed device in figure 2: 



When 2 2 



and when I z 



2:j, 



!U<«)+£-0; 



ttz- 



-Vg(z)=0. 



(36) 



(37) 



The temperature in the cell windows is v 8 (z) and 
the temperature in the gas is V(,{z). The power ab- 
sorbed per unit volume in the cell windows is 

(Ji= a*W '\\uf\. The absorption coefficient of the win- 
dows is a s . 

The boundary conditions are that 



V s (±Z-,)=To, 



(38) 



r)Eo(^i, ru r)rdr. (32) where To is the ambient temperature. Continuity of 

temperature and conservation of heat flow across the 
solid-gas interface at z = ±z 2 give, respectively, 
Appendix A contains the evaluations of the several 
integrals which occur in eq (31) and eq (32). Using eq v s (±z 2 ) = v,,{±Z2) (39) 
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and 



Ksdvsldz | 



: Kffdvu/dz 



(40) 



We integrate eqs (36) and (37) twice to obtain the 
general solutions, 



and 



vAz) =- (Qiz*l2K s ) +C s \z\ +Z) S , (41) 



v M (z)=Cg\z\+D g9 (42) 



where the C's and Z)'s are constants of integration. 

The absolute value of z occurs in eqs (41) and (42) 
because v(z) = v(—z) when a s | z 3 — z>> \ <^ 1. Also, 
because dv g (z)ldz | z =o + =C g and dv g (z)/dz \ z =o~ = 
— C g , the constant C f/ must be zero. (If C g were not 
zero, then a heat sink or heat source would exist at 
z=0. This is not the case for the proposed device). 

The boundary condition (38) at z=z% and the bound- 
ary conditions (39) and (40) at z=z 2 yield three in- 
homogeneous equations for the three unknown 
coefficients C s , D s , and D g ; namely, 



-D ! ,+C S Z2 + D S = (Qiz 2 2 l2K s ) , 
C s K s =Qlz 2 . 



(43) 



A solution for the three equations exists if the deter- 
minant of their coefficients does not vanish; i.e., 




—Ks*0. 



(44) 



Using the matrix method to solve the set of equations 
(43), we ohtain the steady-state temperature: 



v x {z) = T + (Q,IK S ) [z,{ (z-,12) -*} - | z | 

{(|z|/2)-z 2 }], 

S |z I ^ Zs, and 



fo 



for 



v g (z) = T +(Qil2Ks)(z s -z 2 ) 2 



(45) 



(46) 



5. Transient Solution — Plate 

The transient solution u(z, t) for times t 2* satisfies 
the following heat diffusion equations: [5] 

— u s (z,t)=- (47 

dz- k s ot 



for z 2 



:Z :{ 



, and 



— 2 u g (z,t) = 



1 dUg (Z, Q 

k(j 



dt 



(48) 



for | z | ^z 2 . 

At 'time £=0, the sum of the transient and steady- 
state temperature at the exterior sides of the cell win- 
dows at |z|=Z3 is To, i.e., 



because 



^,.(±z ;? , t) = 0, 



Lim u(z, t) = 



(49) 



(50) 



by the definition of transient behavior and because 
v s (±z s ) = To. 

General separable solutions to eqs (47) and (48) are 



v-here 



u(z, t) = R(z)T(t) 



R(z)=[A exp (i&)+B exp (-i£z)] 



r(t) = exp (-£ 2 kt). 

The constant of separation is £ and the coefficients 
A and B are constants. The subscripts s and g denote 
the quantities A, B, £, and k respectively for the 
windows and for the gas. 

Equation (49) gives us that 



R s (±z 3 ) = 0. 



(51) 



The continuity of temperature across the interfaces 
at |z| = Z2 requires that 



«,(: 



: «,,(: 



and 



£l=(k,JkMl- 



(52) 



In addition, the conservation of heat flow across the 
window-glass interface at | z \ = z 2 requires that 



K s (dR s /dz)\ z = ±z . z = K !l (dRJdz) 



(53) 



Hence, the boundary conditions (51), (52), and (53) 
yield six homogeneous simultaneous equations for the 
six unknown coefficients Af, B>, A g , B g , Af, andBf, 
where the superscripts < and > are, respectively, for 
z<0 and z > 0. Because Af = Bf, Bf = A>, and 
A g = B g by symmetry, these six equations reduce to 
three equations in three unknowns. A S = A>, B S = B>, 
and A = 2A g . Solutions for this set of equations exist 
if, and only if, the determinant D of their coefficients 
vanishes. The vanishing of the determinant yields the 
eigenvalues f #=£;, where 7= 1, 2, 3, . . . °°. 
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The;th eigenvalue f, is that value of fr for which the We multiply both sides of eqs (59) and (60) by 

determinant vanishes for the jih time, i.e., £/<frn! 



namely, 



detfl(fr) = 0, 



(54) and by 



r 

D=\ — COS (gffZ2) 

\ jjyKg Sin QgZ 2 



exp (ifjsZs) 

exp (i^sZ-i) 
i^sKs exp (igsZ-i) 

exp ( — i^.s-z.-O 
exp ( — ifrz^) 
i^sKs exp ( — igfo 



(55) 



am/ £§ = {k<,lk s )£;j r From eqs (54) and (55), the eigen- 
value £j satisfies the following equation 

(K g lkf) sin (frit) sin (eJ)=(KJkll*) 

cos (^,,2 2 ) cos (frd), (56) 
where fr = £j and d= (z ; <— 22). 

Expressing /4 a and B s in terms of A=Aj, we have 



/?,,(z) = £ ^j cos (frz) for |z| ^ 22 



(57) 



and 



j=i Vb .y ' 



•cos (fj(W+*t-2ft))} 



(58) 



for z 2 ^ \z\ ^ 2 :} , where f J = {k 9 lk s )^^j- 



We determine the coefficients ^j by using the £ = 
relations obtained from eq (49); namely, 



and 



Rs(z) = To-v s (z) for 22 **\z I ^2., (59) 



R ff (z) = r -%W for |z| ^22. (60) 



Equations (57) and (58) contain the eigenfunctions 
for the window-gas system; namely, 



E(j, z) 



cos (frz) for |z| ^22 

cos (frz,) {cos (£(|z|-z 2 )) 

-COS (^.(|2| +22-22,))} 



2 sin* (f^d) 



for 22 



2.,. 



(61) 



K g (kslk g yi 2 E(j,z) when-2, 



K 8 (kglk 8 yi*E{j,z) when |*| <z s , 



2, 



and we integrate the resulting expressions respectively 
over the intervals ^ 2 ^ z 2 and 22 ^2^ z 3 . This 

procedure gives us the following expression for the 
coefficients A j : 



= (To-Dg^giUkg) 1 ' 2 j E(i, z)dz 

+ (QtlKs)Ks(kglks) 112 I** [z{(zl2) 

-z,}-z,{(z.J2)-z 2 }]E(i, z)dz 
where 

A(fr, Zi)=Kg{kslkgyi 2 \"E(j> z)E(U z)dz 

■Ks(kglk 8 ) 112 jE(jf z)E(U z)dz. (63) 



(62) 



+ , 



Appendix B contains the evaluations of the several 
integrals which occur in eqs (62) and (63). Using eq 
(B4) gives us that A(f;, £,) vanishes unless fj = f*. 
And when £/ = £*, the factor A(£;, £/) = Aj becomes 
from eq (B5) 



A; = 



Kg{k 8 ik g yi 2 



[z 2 +{K 8 kgd cos 2 {£jZ 2 )IKgk 8 sin 2 (fjrf)}]. 



(64) 



Applying eq (Bl) through eq (B3) repeatedly to the 
right-hand side of eq (62) gives us an expression for the 
coefficients A-,; namely, 

Aj= (KglcQil^KskgAj) {(kslkg)*' 2 sin fox,) 

-(K. cob (6*,)/X„8in (£>*)}. (65) 

Because the computer program in [5] evaluates a 
finite number of terms in such expressions as eq (57) 
and eq (58), the coefficients A; are normalized for 
1 ^/^JV by using the relation (60) evaluated at 2 = 22. 
Replacing Aj with (i\Aj (N) in relation (60), one obtains 
an expression for the normalization factor a\ , 



a.v = — ■ 



(Qid 2 l2K 8 ) 



2 Aj(N) cos (fr*,) 

j = i 



(66) 
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6. Cyclic Heating and Cooling 

The Green's function method offers a convenient 
way for computing the response of the proposed baro- 
thermal devices when the beam of radiation is turned 
on at £ = and is modulated at angular frequency co for 
all times t ^ 0. We express the power absorbed per unit 
volume Qi(t) by the expression, 



Qi(t) = 7](t)iq(ri — r)Qi(l + m cos wi) 



(67) 



where \m\ ^ 1 is the modulation factor for cyclic 
heating and cooling and where co is less than the 
characteristic frequency of the system. The step 
function r\ (x)' is zero for x ^ and 1 for x > 0. 

6.1. Cylinder 

The temperature T(r {) , to) = u(ro, to) +i;(ro) — To, 
satisfies the heat diffusion for the cylinder and the gas; 
namely, for times to ^ 0, 

VgTVo, t^-k-^dTinu to)ldto)=-p(r , to); (68) 

where p(r () , to) = r)(ri — ro) (QilK 8 ) (1 + m cos cot) for 
r,^r»; where Vg becomes {(d 2 /dr 2 ) + ro^d/dro) 
because the temperature does not depend upon angle 
do and coordinate Zo\ and where k is subscripted with 
s when ^ r () ^ r 8 and with g when r s ^ r<> ^ r\. We 
subtract the ambient temperature To from the tem- 
perature for later convenience in using the boundary 
condition T(r,. t) = 0. 

The Green's function G(r, t\r , t ) for the system 
satisfies the equation 



where k and K are subscripted with s when < r ^ r s 
and with g when r s - ^ r ^ n. 

We shall now express the solutions of the inhomo- 
geneous diffusion eq (68) for the solid and gas in terms 
of the Green's function G. We first multiply eq (68) by 
G(r, t\ro, to) and eq (72) by T(r {) , to) and then subtract 
the first product-equation from the second product- 
equation. Keeping in mind that the thermal factor for 
the solid (K s /k 8 ) differs from the thermal factor for the 
gas (Kylkfj), we integrate the above result over the 
time domain ^ to ^ t— * and over the volume of the 
solid-gas system. 

IV(. . .)=r*(. • -)(K s lks)rodro 



+ 



/'■(...)« 



lky)rodro, 



where t + = lim t + e. 



The differential form of Green's theorem for our 
special case, 

r (JVlG-GVlT) = j-(Tr ^) 
dr,> \ dro/ 



d ( r 3T 



(73) 



(69) 



enables us to express the temperature in terms of the 
Green's function. The final result becomes 

7 ] (r)7 1 (r s -r)T s (nt)^ V (r-r s )i 1 (r i -r)T !l (r,t) 

[t+ % g fr 8 

= dto~r~ p(r , t {) )G{r, t\r {) , t {) )r {) dr {) (74) 
Jo a\v Jo 

-f dto-r 1 p(r () , to)G(r, t\r {) , to)rodr , 
Jo ky Jr s 

and with g when r s ^ r ^ n. We interpret G to be the 

temperature which occurs at r and t when an impulse where we have used the following initial conditions 
source of heat is introduced at r and to. and boundary conditions: 

The Green's function satisfies a causality condition: 



V 2 G(r, t\r , to) ~ k~ l dG(r, t\r , t )ldt 

8(r-ro)8(t-to) 
r(Klk) 



where k and K are subscripted with s when 0?r? r x 



G(r,t\r<>, t ,) = 



if t< t». 



(70) 



In addition, it satisfies a reciprocity condition under 
time reversal because the causality condition (70) 
requires a time sequence 



G(r, t\r , to) = G(ro, — to\r, — t) , 



(71) 



7,(r, 0) =0, T„(r, 0) = 0, C(r, t\r , t+)=0 

T s (r 8 , t) = T g (r 8 , t), T g (n, t) = 0, 

k~ x G{r, t\r~, to) = k~ l G(r, t\r+,t ), 

G(r, t\n, *o)=0, 
Kgd TJd r | ro = r 8 = kgd TJd r 1 r„= r s , and 



where the function G(r , — to\r, —t) gives the tempera- 
ture at ro and -to due to a source impulse at r and K^dG{r, t\n, t )ldro\r Q =r 8 
— t. Because —to> — t, the time sequence is properly 

= K {l k- ] dG(r, t\r , to)/dr | r ,=rJ. 



ordered. Hence, from eq (71) we have that 
V?)G(r, t\r , t ) + k~ l dG(r, t\r , t )ldt 

8(r-ro)8(t-to) 

ro(Klk) 



We now proceed to construct the Green's function 
(72) C( r > */ r <>, to). Referring to eqs (30), (32), and (33), we 
have that the spatial eigenfunctions of the solid-gas 
system are 
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Wj(r)-- 

where 



(Kslk 8 Ajyi 2 Eo(gj , n , r) for *£ r ^ r, (75) 
(KJk.Aj) [ i 2 Eo(£] , r« , for r,. ^ r ^ n , 



^ tv * ^ / i dT{zo, to) 

— 1 (zo, to) — k~ x : = — p(zo, to), 

azr, dt 



(81) 



f 



rdrW jl (r)W h (r) = 8 h ,j 2 ; 



where p(z , £<>) = tj( |z| — z 2 ) 7?(z : < — |z|) (Qi/K s ) (1 
+ /n cos wO and where £ is subscripted with s when 
z 2 ^ \z\ ^ z :} and with g when z 2 ^ |z| ^ 0. 

The Green's function £(z, j|z<>, to) for the gas cell 
satisfies the equation, 



2 lTj(rWj(ro) = {v(r) 7,(r,-r) (KJk x ) 



d'Gldz--k-'dGldt-- 



8(z-Zo)8(<-n 
K(k'lk)"- 



+ V (r-r s ) V (r i -r)(K,,lk !l )}r-<8(r-r l) ) 
and V 2 r j (r)+^jr j (r) = 0. 



(82) 



Hence, we may expand the Green's function in terms 
of the above eigenfunctions Wj(r); namely, 



where, when 0^ \z\ ^z 2 , k=k u , k' = k 8 , and K=K g ; 
and where, when z 2 ^ \z\ ^ z :{ , k= k s , k' = k ih and 
K = K S . This Green's function for the gas cell obeys 
causality and conditions similar respectively to eq 
(70) and eq (71). It also satisfies the adjoint equation 
similar to eq (72). 

Expressing the temperature in terms of the Green's 
function involves the same steps as those contained in 
part A of this section. The thermal factors are 

Substitution of eq (76) into eq (69) yields the first K s (k,,lk s ) 112 for the windows and K (l {kjk u ) v - for the 

order differential equation for the coefficients C/ (£, to), gas. The final result is, 



G(r,*|r ,*o)= 2 Cj(t,t Wj(r)Wj(ro). 



(76) 



k-*-Cj(t,u ) ) + eC j (t,t») = 8(t-to). (77) T g (z, = 7*0+ VAj cos (£jz)f(a>, £j,k y , m, t), (83) 

fit J ^" 



dt 



j = i 



The solution of eq (77) subject to the condition that where Aj is given by eq (65), /(co, £,-, kg, m, t) is given 

by eq (80), and where we have added the ambient 
temperature To. 



Cj(t, t {) ) =0tot< to is 

Cj (t,to) = 7)(t— to ) CXp {- (*kg (t -to)}. 



(78) 



The temperature in the gas {r 8 ^r^ r e ) for the 
cyclic heating and cooling described by eq (68) becomes 
according to eq (74) 

T g (r, t) = T + J AjFoUj, r u r)f(a>, &, k 9 , m, t) 



; = i 



(79) 



7. Conclusions 

The mathematical expressions contained in eqs (33), 
(34), (35) and (79) for cylinders and in eqs (64), (65), (67) 
and (83) for plates form the basis of the computer 
programs used in references [3] through [6], The 
authors of these references predict by computer pro- 
grams the performance of the proposed devices for 
measuring absorption coefficients as functions of their 



u u jj j ,u u • + , t operating parameters, 

where we have added the ambient temperature io, 

where 



f(a>, £j, k 9 , m, t) = [{exp(-^A,0-l} ( 80 ) 

-f ra{l + (col^ 2 k u ) 2 } " ' { exp (— ijlkgt) — cos cot 



8. Appendix A. Integrals Occurring in 
Equation (31) and Equation (32) 



where Aj is given by eq (34). 



6.2. Plate 



If Z?„(£r) is anyone of the functions J tl (£r), F„(£r), 
or F„(£;, n, r), then from references [9], [10], and [11], 
- ((ol^kg) sin cot}] , we hav e that 

/ Bo(£r)rdr= (r/f )J,(fr) + constant, (Al) 
f Bo(r)r»dr= (rl£)[{r*- (2/0 2 }B>({r) 

+ (2rl£)Bo(£r)] + constant, (A2) 



The temperature, T(z , to) = u(zo, t)+v(zo)-T , f ln (r)Bo(£r)rdr= (r/f) In (r)£,(£r) 
satisfies the heat diffusion equation for the gas cell; 
namely, for times to ^0, + (l/£ 2 )#<>(£r) -f constant. (A3) 
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If Z m ({ir) = aJ m ({ir) + bY m (€ir) and X m (fyr) 

= cj m {^ir) -\-dY m (^ 2 r) with a, b, c, and d independent 
of m and r, then 

{&-&)$ Zm^r)Xm{frr)rdr 

= €irZm(€ir)X m -ifar)-girZm-i(i;ir)Xm(frr). (A4) 

9. Appendix B. Integrals Occurring in 
Equation (62) and Equation (63) 

We tabulate in this Appendix the several integrals 
which appear in eq (62) and eq (63): 



/ cos {£(z+b)}dz = (Hi) sin {£ (z+6)} 
J z cos {%(z+b)}dz=(\lZ) 2 cos {£(z+b)} 

+ (-z sin {£(z+6)} 
JV cos {£(z+&)}dz=(l/£) : *[2£z COS 

{£(*+&)}+ (fV-2) rintt(2+6)}]. 
When f 1 7^ ^2, we have that 

/ cos Ui(z+b)} cos{f>(z + c)} 



(Bl) 



(B2) 



(B3) 



8in{(g,-f,)U+6)} 
2(fi-6) 



: COS {^2(^~ C)} 

sin{(g,+&)(z-t-6)} 



+ ■ 



2(£i + £i 
+ sin {£2(6 -c)} 



c»s (£. -&) (z+6) 
2(^,-6) 



cos{(g, + g*)(z+6)} ] 

2(^ + 6) J' (B4> 

and when f j = £ 2 = £ , we have that 

Jcos{£(z+6)} cos{£(z+c)} 

_ cos {g(6 — c)} 
2£ 

[£(z+6) + sin {£(z+b)} cos {£(z+6)}] 

+ sin{^-c)} s . n2{ ^ +6)}j (B6) 
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